Adiabatic quantum computation: Enthusiast and Sceptic's perspectives 



Zhenwei Cao and Alexander Elgart 
Virginia Tech 
(f^Dated: May 3, 2010) 

Enthusiast's perspective: We analyze the eftectiveness of AQC for a small rank problem Hamil- 
tonian Hp with the arbitrary initial Hamiltonian Hj. We prove that for the generic Hi the running 
time cannot be smaller than 0(y/N), where A'^ is a dimension of the Hilbert space. We also construct 
an explicit Hi for which the running time is indeed 0{^/N). Our algorithm can be used to solve 
the unstructured search problem with the unknown number of marked items. 

Sceptic's perspective: We show that for a robust device, the running time for such Hp cannot be 
much smaller than 0(A^/lnA''). 

PACS numbers: 03.67.Lx, 03.67.Ac, 03.65.Aa, 03.65.-w, 02.10.Yn 



Overture. — Adiabatic quantum computation (AQC) 
(e.g. [Tj) is a Hamiltonian-based model of quantum 
computation. The idea behind AQC is that finding the 
ground state of a problem Hamiltonian Hp solves inter- 
esting computational problems. In the abstract setting, 
let H] , Hp be a pair of hermitian N x N matrices, with 
A^ ^ 1. Consider the interpolating Hamiltonian H{s) of 
the form 

His):^{l-f{s))Hi + fis)Hp, (1) 

where / is a monotone function on [0, 1] satisfying /(O) = 
0, /(I) = 1. We will denote by Pj (respectively Pp) 
the spectral projection on the ground state energy Ej 
(Ep) of the matrix Hj (Hp). We prepare the initial 
state of the system ip{0) in the {a-priori known) ground 
state tpi € Range Pj of the Hamiltonian Hj p] , and let 
the system evolve according to the (scaled) Schrodinger 
equation: 

tMs) = tH{s)Ms), Vr(0) = (2) 

The adiabatic theorem of quantum mechanics ensures 
that under certain conditions the evolution ipri^) of the 
initial state stays close to the Range Pp. For AQC 
to be a potent quantum algorithm, the running [i.e. 
physical) time r in Eq. ([2| must be much smaller than 
N. Although AQC attracted a considerable interest in 
physics and computer science communities, the quantita- 
tive characterization of the speed up in its use remains at 
large unknown. The core issue here is related to the ex- 
treme sensitivity of the adiabatic behavior to the spectral 
structure of the operator H{s). Specifically, the devia- 
tions may become large when the gap g{s) between the 
ground state of H{s) and the rest of its spectrum is small 
in the vicinity of some instant s G [0, 1]. 

The traditional approach to the problem so far was 
to estimate this minimal gap [5]. Putting a few rare 
exceptions aside {e.g. [1]), it is usually a very hard task. 
This explains why, generally speaking, not much light was 
shed on the effectiveness of AQC. Let us note that the 
estimates of the running time involving the gap g{s) alone 



provide only the upper bound on the optimal running 
time r. In reality r can be much smaller. 

In this paper, we discuss the reliable upper and lower 
bounds on the optimal value of r, circumventing the es- 
timates on the size of the gap. Our method is applicable 
for a particular class of problem Hamiltonians, satisfying 
the following hypothesis. 

Assumption 1. The problem Hamiltonian is of the 
small rank: Rank{Hp) :— m ^ N . 

Even in this narrower context, there is no unequivocal 
riposte to whether AQC is indeed efficient, as we shall see. 
As often happens in theoretical deliberations, the answer 
depends, to some extent, on the degree of your zeal. To 
keep the discussion balanced, we present two different 
perspectives: The first one is on the optimistic side while 
the second one is rather pessimistic in its nature. To this 
end we set a stage for two close acquaintances, Messrs. 
Enthusiast and Sceptic, and let the wise Reader judge 
who of them is closer to the mark. 

Let us note that for AQC to work, it suffices to ensure 
that V'r (1) has a non trivial overlap with the range of Hp, 
which we will encode in the requirement ||Pf'0t-(1)|| > 
1/5, [S]. Another issue that usually arouses certain de- 
gree of confusion, which we want to avoid, is a normal- 
ization of H(s). To that end, we will use the calibration 
\\Hi\\ = \\Hp\\ = 1. One should bear this convention in 
mind when performing comparison with other results. 

The rest of the paper is organized as follows: We 
first present the discussion from Enthusiast and Sceptic's 
points of view, indicating briefly the intuition behind the 
corresponding assertions. We then give proofs of Theo- 
rems |3] and [4] (the rest of the proofs can be found in [5]). 
Now we pass the baton to Mr. Enthusiast. 

Enthusiast's perspective. — To formulate the result, let 
me introduce a set of the related parameters. First, I 
want to quantify the overlap between the initial state "0/ 
and the problem Hamiltonian. Namely, let = jjiJ^V'/llj 
let 62 = \\Pptpi\\, and let ^3 — \\Qpipi\\, where Qp is a 
projection onto Range Hp. Note that for a generic Hj 
all (5's are small, with Si and ^3 being 0{yJm/N), while 



2 



52 = 0{yjm' /N). Here m' is a dimension of Range Pp- 
Second, I want to distinguish between Ep and the rest of 
the spectrum of Hp^ which I will assume henceforth are 
separated by the gap gp. Finally, since I don't want to as- 
sume that Hp is sign definite and given that \\Hp\\ = 1 by 
convention, the energy Ep will show up in the estimates. 
The prototypical example covered by our results is the 
generalized unstructured search (GUS) problem, which 
can be cast in the following form: Suppose Hp is diago- 
nal with the unknown number of entries equal to —1 and 
the rest of the entries equal to zero (so that Hp — 1—Pp). 
Pick Hi = -\iPi){'iPi\ with V/ = N-^/^{l, . . . , 1). Then 
the corresponding parameters are 63 — 62 — y^m/N, 
Ep = — 1, and gp = I. 

The pair of results below, coupled together, gives fairly 
tight lower and respectively upper bounds on the optimal 
running time in AQC. 

Theorem 1. Consider the interpolating family Eq. (1) 
with an arbitrary f. Then the running time t_ in Eq. (2) 
for which \\Q ptpr^ > 1/5 satisfies 
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for §2 < 1/5 . 



(3) 



The quantitative measure of how much ■(/'t(s) devi- 
ates from tpi is encoded in the size of the commutator 
[Pi,H{s)]. Hence one expects to see the deviation from 
ipi over the time r such that r • || [P/, _ff (s)]|| — 0(1). 
Since Pj commutes with Hj while < 2Si, we 

get \\[Pi,H{s)]\\ < 2Si for all s and the bound in Eq. ^ 
follows up to a constant. 

Let me note that the similar, albeit less sharp (with 
the wrong dependence on m) lower bound was recently 
established in [7j- 

Theorem 2. Suppose S^/gp ~ 0(l/lnA^), Then there 
exists an explicit rank one Hj and an explicit function f 
such that ||PFV'r+(l)|| > 1/5 for 



C{1 - Ep) 
lEplSo 



(4) 



for any C E [1/3,2/3] 



For N ^ m the requirement on 63/ gp is typically sat- 
isfied. Note also that t^/t^ ~ \Jra' jm. This is not 
particularly surprising, as in Theorem [T] the aim was to 
ensure that iI)t_ (1) has an overlap with the range of Hp^ 
whereas in Theorem [2] we want V''r+(1) to overlap with 
Range Pp. 

The choices in the theorem are: Hj = and 
a (non adiabatic) parametrization /(s) is given by 



fis) 



0, 
a 

1, 



s = 

s e (0,1) 

s = 1 



That means we move diabatically (instantly) to the given 
point of the path, stay there for the time t+, and then 
move quickly again to the end of the path. Such / is in 
fact optimal for the Grover's problem. 

The intuition behind this assertion is as follows: With 
the above choice for /(s) 



^r(l) = e 



_ -iaT+-{EFPi +Hf] 



Note now that the ground state energy of EpPj matches 
that of Hp and differs from the energies of its excited 
states. Let X be a subspace spanned by vectors in the 
ranges of Pj and Pp, and let X-^ be its orthogonal com- 
plement (so that X X-^ is the whole Hilbert space). As 
usual in adiabatic setting, the transitions between X and 
X-^ are suppressed due to fast oscillations caused by the 
energy differential. Therefore the initial state ipi slowly 
precesses in the X subspace, and by choosing the right 
value for r+ one can find the evolved state sufficiently 
close to Range Pp. The argument identical to the one in 
Theorem [T] shows that the running time is roughly 

1- Ep 1 _ l-Ep 1 

\Ep\ ■ \\[Pi,Pp]\\ ^ \Ep\ '2h' 

Since the precession is very slow, t+ is fairly robust. 

This assertion can be seen as an extension of the clas- 
sical result of Farhi-Gutmann [8 on the Grover's search 
problem. For GUS the parallel result was established for 
the quantum circuit model (QCM) in ^ . 

Theorem [2] uses values of Ep and 62 = \\PpiJi\\ — 
0{1/VN) as the input. In many important applications 
(such as GUS) the value of 82 is unknown. To this end, 
we prove the following assertion. 

Theorem 3. Suppose that the value of Ep is known. 
Then there is a Hamiltonian - based algorithm that de- 
termines 62 with 1/iV^ accuracy and requires 0((lnA^)^) 
of the running time. 

Note that the combined running time in Theorems 
[2] and [3] remains 0{\fN). The algorithm used in the 
proof is inspired by the mean ergodic theorem and makes 
use of the fact that the survival probability cp{t) ~ 
('0/|e**^*'-0/) is directly measurable in AQC framework. 
For GUS this problem is known as quantum counting and 
was analyzed in QCM framework in [3]. 

Enthusiast's summary. — Theorem [l] tells us that for 
a generic Hj the running time cannot be smaller than 
0{^/N) . Theorems [2] and [3] construct the explicit Hi 
and the parametrization /(s) so that t = 0{\/N). I 
have assumed that the ground state energy Ep of Hp is 
known with the 1/iV accuracy. 

Sceptic 's perspective. — Let me first point out two short- 
comings of the method which is usually employed in es- 
timation of the running time of AQC {e.g. |3] for the 
Grover's problem and [10] ). The technique hinges on 
a choice of a parametrization / such that /(s) is small 
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whenever the instantaneous spectral gap g{s) is small 
[TT] , To construct such /, one need to know the values 
Sj for which g{sj) = 0{N^^/^) with high precision. Such 
analysis requires the detailed information about the spec- 
tral structure of Hp. The similar issue is present (albeit 
to a lesser extent) in the Enthusiast's approach, as one 
still needs to know Ep. Even if this technical hurdle can 
be overcome, the extreme susceptibility of V'r(l) to the 
parametrization / poses a radical problem in practical 
implementation. Indeed, it is presumably extremely dif- 
ficult to enforce / = for a long stretch of the physical 
time, as the realistic computing device inevitably fluctu- 
ates. So in the robust setting one can assume that for 
any given moment sq the value /(sq) is greater than some 
small but fixed k. We can then as well consider the func- 
tions / in the robust setting that satisfy f{s) > k for all 
values of s. 

To understand how the robust system evolves, let me 
consider the following semi-empiric argument, substanti- 
ated in Theorem [4] below. One can show that for a 
finite rank matrix Hp the minimum value g of the gap 
g{s) between the ground state energy E(s) of H{s) and 
the rest of the spectrum of H{s) is 0(^3). Let g = g{so), 
then one can introduce two different time scales: ri and 
T2. We set Ti = A~-^, where A is a gap between the two 
smallest eigenvalues of H{so) and the rest of its spec- 
trum. The scale T2 is associated with a two level system 
corresponding to the restriction of the Hilbert space to 
the spectral subspace of these two eigenvalues. Typically, 
Ti T2, and ipT stays close to the ground state provided 
r » r2. If Ti ^ T <C r2, then ipr will still stay close 
to the range of the above spectral subspace. However, it 
will behave as if the avoided level crossing is a true level 
crossing, with evolution following the first excited state 
rather than the ground state (see Figure [l]) . To estimate 
T2 let me consider a two level system of the form 



Kis) 



l-f{s) Sf{s) 
Sf{s) f{s) 



(5) 



Assuming that / is differentiable, the Landau-Zener for- 
mula shows that for K(s) the corresponding time scale 
is T2 ~ /(so) . The value of the minimal gap here is 
equal to S\/l + (5^ w S. Since for H{s) the value of g is 
roughly 83, we see that tAt stays close to the ground state 
of H{s) only if 

T > f{so)-{S3)-' = fiso)-OiN). 

Hence in the robust setting, where / cannot be too small 
at any given instant, there is no significant speed up in 
using AQC. The result below reaffirms this argument 
for the case of the initial Hamiltonian Hj of the small 
rank. In what follows, Qj will denote the projection onto 
Range Hj, with S := ||Q/QfI1- 

Theorem 4 (Robust lower bound on the running time). 
Suppose that f in Eq. ([T]) is differentiable and satis- 



o-mt)) 




FIG. 1: Illustration for Theorem [i] iV = 10*. Top: A 
pair of relevant eigenvalues of H{t) ~ Hi -f 2tHF as a 
function of t, for Hp = — |ei)(ei| — ||e2){e2| and Hi = 
|<^i){<^ih|02){02|. Hereei = (1, . . . , 0), 62 = (0, 1, 0, . . . , 0), 
(t>i = l/%/iV(l, ■•■,!), and <t>2 = l/%/iV(l, . . . , 1, -1, . . . , -1), 
with exactly half of —1 so that {(j!)i|<^2) = 0. Bottom: The 
magnified region of the avoided crossings in the upper panel. 



fies f{s) > K > for s G [0,1]. Then, if t < Tr 



0\ 



1252 \nS>' 



we have 



(6) 



Hence the running time t for which \\QpijjT-{l-)\\ > 1/5 
cannot be smaller than Tj.. 

Sceptic's summary. — Theorem [4] tells us that for a 
generic Hi of the small rank the robust running time 
cannot be smaller than 0(iV/lnA^). Hence AQC is not 
really effective for the problem Hamiltonians that satisfy 
Assumption [T] 

Proof of Theorem The starting point is a pair of iden- 
tities, ini: 



p^ sin(fca;) 



pP(cosw-l) 



fc=i 



k\ 



sin (p sin w) 



k=l 



p cos(fcw) 
~fc! 



""^^ cos(psinw) . (7) 



In particular, if 1 — cos a; > A, each term in Eq. ([7| is 
bounded by e~^^ and therefore is smaller than pro- 
vided p — 21nA^/A. On the other hand, the remainders 
to the partial sums (up to k = L) in Eq. ^ are 0{p^/L\) 
provided the latter quantity is small. Combining these 
observations, we get 



^ kl 

fe=i 



l + Oil/N^), w = 
0{1/N^), 1-cosw > A ' 



for p ^2 \nN/A and A < 1. Hence 
^"''E^^^^l^'*^''""''"'^^) = {S2f + 0{1/N'), (8) 
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for p = 21niV/(l — cos gp). The total running time is 
ELii-0((ln7V)2). □ 

Proof of Theorem^ For a solution ipr{s) of let 
0,(s):= e»''(^)-^^V,(s), r(l-/M)dr. (9) 



As a result, we obtain a bound 



Then one can readily check that (pris) satisfies IVP 

iMs) = tH{s)Ms), MO) = (10) 

where H{s) = (1 - /(s)) {Hi - Ej) + f{s) Hp- Clearly 
= \{4'i\M\- Let 

= - ((1 - f{s))Ei + ez))-i , 

and let 0(s) = ipi — f{s)HpB{s) -07, where e is a small pa- 
rameter to be chosen later. Omitting the s dependence, 
we have 

H4>^~f{l- fWiHpBiji + iefHpB^i. (If) 

That means that away from the m values of s for which 
B[s) has zero eigenvalue, ||i?0|| is very small, since 
WHiHpiljiW < S^. Note now that 



m)\M\)) = m)\'\>rm 



But (0(O)|</>^(O)) = 1 and 



1 d_ 
ds 



(12) 



I(0(l)l0r(l)>| = 



(V;/|0.(1)) - (V/ 



:|'^r(l)> 



< |(V/|'/'r(l))| + IIQfV'/II < l(^/|0r(l)>| + ^ 



Substitution into Eq. ( f 2 ) gives 

1 - |(0/|0r(l))| < /' 

Hence Eq. ([6]) will follow if 
''^((/.(s)|(/..(s))ds 
We have 



(0(s)|0,(s))ds 



+ <5. 



< 1 



2^6 



2^. (13) 



(</.(s)|<^,(s)) = (0(s)|0.(s)) - zr((/.(s)|i7(s)|(/.,(s)) 



We bound the first term on the right hand side by 
and the second one by Tljif^T-l!- A straightforward com- 
putation (using Eq. (11 1 for the second term) shows that 



A, 



where A,(s) := dist{f{s)a{Hp) , {1~ f{s))Ei+ei). Here 
dist(S', z) is an Euclidean distance from the set S to the 
point z in C, and a{H) stands for the spectrum of H. 



< {f + rS + Te)-^ + 2^^ 
" Ae {A,f 



Integrating both sides over s and using the bounds 



^ fds , fds 2m 
< -2mlne: / — ^ < 



A,(s) 
1 ds 



{AMr 



Inf 



^ , < -2m- 

Ae[s) K 



we can estimate 



' d_ 
ds 



{^{s)\Ms))ds 



< 2mS -Ine 1 



t6 re 



Hence the required bound in Eq. ( 13 ) follows with the 



choice e = 10 "^mS, provided r < — where C is a 



constant. 



□ 
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